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THE  SHOCK  HUGONIOT  OF  MINERAL  OIL 


ABSTRACT 

The  shock  Hugoniot  of  heavy  mineral  oil  (density  =  0.87g/cc)  has 
been  determined  by  plane  shock  wave  experiments.  A  charged  capacitor 
technique  was  used  to  measure  shock  velocities  through  specimens  of 
mineral  oil.  Shock  velocities  through  polymethyl  methacrylate  reference 
specimens  were  measured  by  a  shock-induced  polarization  technique,  and 
impedance  calculations  were  performed  to  establish  the  Hugoniot  equation 
of  state  in  mineral  oil.  Within  the  pressure  range  from  15  to  150  kilo- 
bars,  the  Hugoniot  curve  of  heavy  mineral  oil  was  found  to  be  represented 
by  the  linear  relationship  U=2.19  +  1.52u,  where  U  is  shock  velocity, 
u  is  particle  velocity,  and  velocity  units  are  mm/ysec. 
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I.  INTRODUCTION 


Mineral  oil  has  been  used  in  shock,  wave  experiments  to  eliminate 

1* 

electrically  noisy  air  shock.  It  is  particularly  well  suited  to  this 

application  because  it  produces  a  very  small  shock-induced  electrical 
2 

signal,  and  therefore  does  not  contribute  to  other  electrical  signals. 

The  Ballistic  Research  Laboratories  have  also  used  heavy  mineral  oil 
3 

(UMO)  to  improve  impedance  matching  at  interfaces,  and  as  a  shock  trans¬ 
mitting  medium.  In  the  latter  application,  the  Uugoniot  of  UMO  was 
needed  to  establish  shock  pressures. 

The  Hugoniot  of  UMO  was  determined  using  the  impedance  match  tech- 
4 

nique.  In  this  technique,  specimens  of  UMO  and  polymethyl  methacrylate 
(PMMA)^  were  placed  on  a  metal  plate  through  which  an  explosively  pro¬ 
duced  plane  shock  wave  was  transmitted.  Transit  times  of  shock  waves 
through  the  PMMA  reference  specimens  were  measured  from  the  shock-induced 
polarization  signal*^  and  used  to  calculate  the  average  shock  velocities. 
The  polarization  signal  from  UMO  is  very  small,  so  a  charged-capacitor 
technique^  was  used  to  measure  the  shock-wave  transit  times  from  which 
average  shock  velocities  were  calculated.  Uugoniot  points  for  HMO  were 
then  obtained  by  impedance  calculations.  This  report  describes  the 
experimental  details  and  presents  the  Hugoniot  data  obtained  for  HMO. 

II.  EXPERIMENTAL  TECHNIQUES 

A.  Shock  Pressures 

Plane  shock  waves  were  generated  by  the  explosive  train  shown  sche¬ 
matically  in  Figure  1.  This  explosive  train  consisted  of  a  10-cm  dia¬ 
meter  composition  B-TNT  plane  wave  lens  and  a  10-cm  diameter  by  2.5-cm 
thick  base  charge  of  TNT,  composition  B,  or  75/25  Octol.  The  plane 
shock  wave  from  the  explosive  was  transmitted  to  the  test  specimens 
through  a  buffer  plate.  The  higher  shock  pressures  were  produced  by 
using  different  base  charges  on  0.6-cm  thick  buffer  plates  of  AZ31B 
magnesium,  2024  aluminum,  or  A.S.T.M.  B16  brass.  The  low  shock  pressures 

*Rei$eAence2>  <v ie  &ound.  on  page  25 
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Figure  1.  Experimental  Arrangement  for  Simultaneous  Measurement  of  Shock 
Velocities  in  Mineral  Oil  and  Polymethyl  Methacrylate 
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were  produced  by  using  a  TNT  base  charge  on  a  laminated  buffer  which  re¬ 
duced  the  pressures  by  impedance  mismatch  between  laminations.  The  lam¬ 
inated  buffer  consisted  of  three  0.6cm  thick  laminations,  brass-x-brass, 
where  x  was  low  density  polyethylene  or  AZ31B  magnesium. 

These  systems  produce  shock  waves  which  are  plane  to  0.1  micro¬ 
seconds  over  diameters  of  two  to  three  and  one- half  inches.  Each  experi¬ 
ment  is  designed  to  use  only  the  plane  portion  of  the  wave.  The  plane 
wave  is  delivered  normally  incident  to  the  buffer-sample  interface.  In 
practice  some  of  the  waves  will  arrive  with  measurable  obliquity.  This 
is  assessed  and  compensated  for  by  measurements  of  signal  risetime,  as 
described  in  Section  III. 

B.  Experimental  Arrangement  for  Mineral  Oil  Measurements 

Shock  velocity  through  mineral  oil  was  measured  by  a  charged  capac¬ 
itor  technique.  The  test  capacitor  shown  in  Figure  2  consisted  of  a 
1.27-cm  diameter  brass  second  electrode  suspended  0.10-cm  from  the  buffer 
plate  which  served  as  the  first  electrode.  The  HMO  sample  was  introduced 
into  the  interelectrode  gap  and  was  retained  by  capillary  action.  The 
test  capacitor  was  charged  to  a  potential  of  475  volts.  A  0.01  micro¬ 
farad  capacitor,  in  parallel  with  the  test  capacitor,  held  the  applied 
voltage  essentially  constant  during  an  experiment.  The  test  capacitor 
was  connected  in  series  with  a  coaxial  cable  terminated  by  its  character- 

g 

istic  impedance  at  the  input  of  a  fast  rise  oscilloscope.  Under  static 
conditions,  no  current  flowed  in  the  circuit  after  the  capacitor  was 
initially  charged.  When  the  explosive  charge  was  detonated,  a  plane 
shock  wave  passed  through  the  buffer  plate  and  entered  the  HMO,  com¬ 
pressing  it  and  increasing  the  dielectric  constant.  This  increased  the 
capacitance  of  the  test  capacitor  and  caused  a  charging  current  to  flow 
in  the  circuit.  The  resulting  voltage  drop  across  the  terminating 
resistor  was  recorded  with  an  oscilloscope. 

As  shown  in  Reference  7,  the  signal  profile  may  be  calculated  using 
the  equation 
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(1) 


Ve0kQA[U-(k0/k)(U-u)] 

{  XQ- [U- (kQ/k) (U-u) ] t } 2 

where  V  is  the  applied  voltage,  e  is  the  permittivity  of  free  space 

-12  u  ' 

(8.854  x  10  F/m) ,  A  is  the  area,  k^  is  the  relative  dielectric  constant 

at  zero  pressure,  k  is  the  relative  dielectric  constant  of  the  shock  com¬ 
pressed  dielectric,  xQ  is  the  initial  thickness,  and  t  is  the  time  mea¬ 
sured  from  the  entry  of  the  shock  front  into  the  dielectric.  The  dielec¬ 
tric  constant  under  shock  compression  may  be  roughly  estimated  by  the 
Clausius-Mosotti  equation,  assuming  constant  polarizability.  The  Clausius- 
Mosotti  equation  was  used  instead  of  the  Drude  equation,  which  Hauver 
(7)  has  shown  to  be  more  appropriate  for  polyethylene,  because  it  pre¬ 
dicts  a  larger  change  in  dielectric  constant  and  therefore  a  larger 
signal.  This  increases  the  safety  factor  in  the  estimate  of  signal  size 
used  to  set  the  oscilloscope  and  reduces  the  chance  of  an  off-screen  sig¬ 
nal.  It  should  be  noted  that  the  calculation  assumes  a  three-electrode 
parallel  plate  capacitor.  The  simple,  unguarded  capacitor  design  used 
in  these  experiments  is  not  suitable  for  quantitative  measurements  of 
the  dielectric  constant  under  shock  compression,  since  it  provides  no  pro¬ 
tection  from  edge  effects.  Nevertheless,  the  calculation  provides  a  rough 
estimate  of  signal  size,  and  the  unguarded  design  does  provide  signals 
which  may  be  readily  interpreted  and  measured  for  shock  transit  time, 
which  is  the  only  value  needed  for  the  Hugoniot  measurement. 

C.  Experimental  Arrangement  for  Polymethyl  Methacrylate  Measurements 

The  shock  velocity  through  the  PMMA^  reference  specimens  was  measured 
by  a  shock-induced  polarization  technique. ^  The  experimental  arrangement 
is  shown  in  Figure  3.  A  PMMA  specimen,  0.72cm  in  diameter  by  0.10cm 
thick,  was  placed  in  contact  with  the  buffer  plate  which  served  as  the 
first  electrode.  The  other  surface  of  the  PMMA  specimen  was  covered  with 
a  spring-loaded  brass  disc  which  served  as  the  second  electrode.  A 
grounded  aluminum  ring  was  used  as  electrical  shielding.  The  electrodes 
were  connected  to  an  oscilloscope  by  a  coaxial  cable  which  was  terminated 
with  its  characteristic  impedance  at  the  oscilloscope  input.  When  the 
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Figure  2.  Mineral  Oil  Capacitor 


Figure  3.  Polymethyl  Methacrylate  Shock-Induced  Polarization  Transducer 
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plane  shock  wave  entered  the  PMMA  specimen  from  the  first  electrode, 
polarization  was  induced  and  a  displacement  current  charged  the  capaci¬ 
tance  in  the  circuit.  The  voltage  drop  across  the  cable  termination  was 
displayed  on  the  oscilloscope.  The  shock  transit  time  was  obtained  from 
the  duration  of  the  polarization  and  was  used  to  calculate  the  average 
shock  velocity. 


III.  SIGNAL  MEASUREMENT  AND  ANALYSIS 

Typical  signals  from  HMO  and  PMMA  are  shown  in  Figure  4.  Despite 
the  difference  in  the  generating  mechanism,  the  two  types  of  signals 
appear  similar  and  were  measured  the  same  way. 

A  sharp  initial  voltage  rise  occurred  at  time  t^  when  the  shock 
wave  entered  the  specimen.  The  risetime  t^-tg  resulted  from  tilt  and/or 
curvature  of  the  incident  shock  front.  After  time  t^,  the  signal  profile 
increased  slowly  until  time  t^,  when  the  shock  wave  arrived  at  the  sec¬ 
ond  electrode.  The  apparent  transit  time  was  correctet*  f°r  rise¬ 

time,  and  was  used  with  the  measured  thickness  to  calculate  the  average 
shock  velocity  through  the  specimen.  Since  the  shock  wave  undergoes 
attenuation  as  it  passes  through  the  samples,  the  calculated  velocity  will 
be  slightly  lower  than  that  at  the  buffer-sample  interface.  The  HMO  and 
PMMA  samples  are  the  same  thickness  and  the  shock  waves  in  both  samples 
will  be  attenuated  by  comparable  amounts.  The  net  effect  will  be  to 

shift  the  measured  point  a  short  distance  down  the  U  -u  curve.  The 

s  p 

attenuation  in  PMMA  has  been  shown  (reference  6)  to  be  about  0.5%/mm. 

The  design  of  the  experiment  therefore  assures  that  any  error  from  this 
source  is  very  small. 

A  risetime  correction  is  necessary  when  the  shock  front  is  obliquely 
incident.  Because  of  pressure  release  at  the  free  lateral  boundary  of 
the  specimen,  times  tg  and  t^  relate  to  shock  paths  that  originate  from 
different  points  along  the  incident  shock  front,  necessitating  a  correc¬ 
tion  to  t2“to*  obliquely  incident  shock  front  is  considered  to  be 
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Figure  4a.  Mineral  Oil  Change  of  Dielectric  Constant  Signal.  Vertical 
scale  0.2  volts/div,  horizontal  scale  0.04  usec/div,  timing 
marks  0.020  usee  intervals. 


Figure  4b.  Polymethyl  Methacrylate  Shock- Induced  Polarization  Signal. 

Vertical  scale  0.2  volts/div,  horizontal  scale  0.04  usec/div 
timing  marks  0.02  usee  intervals. 
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the  more  general  condition  encountered,  so  the  risetime  correction  is 
always  applied  when  the  signal  risetime  exceeds  the  oscilloscope  rise¬ 
time.  The  RC  response  of  the  circuit  is  less  than  the  instrument  rise¬ 
time. 

The  corrected  transit  time  for  a  plane  but  obliquely  incident  shock 
wave  has  been  shown  to  be  (see  Reference  6) . 


T 


(t2-t0)-At(-) 


(2) 


where  T  is  corrected  transit  time,  At  =  t^-tg,  w  is  the  diameter  of  the 
sample  and  r  is  the  curvature  (the  radial  distance  from  the  edge  of  the 
sample  to  the  point  where  the  shock  wave  first  strikes  the  rear  elec¬ 
trode)  . 

The  curvature,  r,  for  PMMA  has  been  measured ^  and  may  be  applied 
directly.  The  curvature  for  HMO  has  not  been  measured  and  was  estimated. 
This  was  done  by  first  calculating  uncorrected  values  for  the  Hugoniot 
equation  in  the  form  U  =a+bu,  where  a  and  b  are  constants.  The  approxi¬ 
mate  sound  velocity,  c,  of  the  lateral  release  wave  at  each  experimental 

9 

pressure  was  then  calculated  using  Jacobs'  approximation 


c  = 


^  CU-bu) 


(3) 


★ 

The  approximate  curvature,  r  ,  was  then  calculated  from  the  geometric 
relationship  of  the  compression  and  lateral  release  waves.  As  shown  in 
reference  6,  the  true  sound  velocity  is  given  by 


c 


2  P0,21l/2 


u[(r-)  +  (r2-)  ] 


where  c  is  the  sound  velocity  and  x^  is  the  initial  thickness  of  the  HMO 
specimen.  This  equation  may  be  written 


c 


u[£  )2  - 

xo 


,U-u, 2. 1/2 
1  U  J  J 
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or,  solving  for  r. 


,,c.2  U-u.2.1^2 

r  ■  vy  -  <-—■ 1 1 

Since  the  value  calculated  for  c  by  Jacobs  equation  is  an  approximation, 

*  * 
c  ,  the  value  calculated  for  the  curvature  r,  will  be  r  ,  the  approximate 

curvature 

*  *  9  „  2  1/2 
*  _  r  z  Il-ii 

(4) 


*  9  2  1/2 
r,c  .2  ,U-U.  , 

'  xo[(; )  -  (— J  1 


The  approximate  curvature  for  the  mineral  oil  was  used  in  Equation  (2) 
to  calculate  corrected  transit  times  and  shock  velocities  for  HMO.  The 
corrections  were  small  [At  (— )<  1% Ct 2~ t0) ]  so  The  approximations  involved 
do  not  greatly  affect  the  results. 

The  Hugoniots  of  PMMA  and  the  buffer  plate  metals  are  well  established 
in  the  pressure  ranges  studied  here.^  The  pressure  and  particle  velocity 
in  the  PMMA  were  calculated  from  the  shock  velocity,  using  the  relation¬ 
ship  U  =  2.695  +  1.538u,  and  the  conservation  relation  P  =  p0  Uu. 

The  method  of  impedance  matching  was  used  to  determine  the  buffer 
plate  conditions.  The  release  adiabat  of  the  buffer  material  must  pass 
through  the  previously  determined  PMMA  point.  This  locates  the  buffer 
pressure  and  particle  velocity  point  and  identifies  the  locus  of  states 
which  can  be  attained  in  a  material  in  contact  with  the  buffer  plate. 

The  reflected  buffer  liugoniot  was  used  as  an  approximation  for  the  buffer 
adiabat.  Since  the  PMMA  and  HMO  Hugoniots  are  in  close  proximity,  this 
is  a  good  approximation.  The  HMO  Hugoniot  point  is  found  at  the  inter¬ 
section  of  the  reflected  buffer  Hugoniot  and  a  line  of  slope  p0  U  (mineral 
oil) ,  where  p0  is  the  initial  density  of  the  HMO. 

The  major  sources  of  error  in  the  experiments  reside  in  the  physical 
measurements  of  the  set-up  components  (estimated  at  ±1%)  .  in  the  measure¬ 
ments  of  the  film  records  of  the  oscilloscope  traces  (estimated  at  ±1%). 
the  effects  of  shock  wave  attenuation  (estimated  to  be  less  than  0.5%) 
and  the  time  error  due  to  shock  wave  obliquity,  which  is  estimated  to  be 
less  than  2%. 
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IV.  RESULTS 


The  experimental  data  are  given  in  Table  I  and  are  shown  graphically 
in  Figure  5.  A  linear  least  squares  fit  of  the  data  gives  the  relation¬ 
ship  U  =  2.18  +  1.53u,  a  =  .076.  A  Hugoniot  data  table  based  upon  this 
equation  is  given  in  Table  II. 
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Figure  5.  Hugoniot  Data  for  Heavy  Mineral  Oil,  Presented  in  the  Shock 
Velocity-Particle  Velocity  Plane. 
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TABLE  I.  Experimental  Data  for  Shock  Velocity  and 
Particle  Velocity  in  Heavy  Mineral  Oil 


xo 

t2"t0 

t2_t0 

U 

u 

mm 

usee 

corrected 

usee 

mm/ usee 

mm/ usee 

0.960 

0.327 

0.326 

2.94 

0.55 

0.94S 

0.270 

0.269 

3.51 

0.86 

0.986 

0.227 

0.226 

4.36 

1.36 

0.940 

0.189 

0.187 

5.03 

1.88 

0.958 

0.168 

0.166 

5.74 

2.32 

0.945 

0.149 

0.148 

6.40 

2.81 
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TABLE  II.  Shock-Wave  Compression  Data  for 
Heavy  Mineral  Oil 


u 

U 

P 

v/vo 

mm/nsec 

mm/nsec 

ki lobars 

.0000 

2.1856 

.00 

1.0000 

.  1000 

2.3377 

2.03 

.9572 

.2000 

2.4897 

4.33 

.9196 

.3000 

2.6417 

6.89 

.8864 

.4000 

2.7937 

9.72 

.8568 

.5000 

2.9457 

12.81 

.8302 

.6000 

3.0977 

16.17 

.8063 

.7000 

3.2497 

19.79 

.7846 

.8000 

3.4018 

23.67 

.7648 

.9000 

3.5538 

27.82 

.7467 

1.0000 

3.7058 

32.24 

.7301 

1.1000 

3.8578 

36.91 

.7148 

1.2000 

4.0098 

41.86 

.7007 

1.3000 

4.1618 

47.07 

.6876 

1.4000 

4.3139 

52.54 

.6754 

1.5000 

4.4659 

58.28 

.6641 

1.6000 

4.6179 

64.28 

.6535 

1.7000 

4.7699 

70.54 

.6436 

1.8000 

4.9219 

77.07 

.6342 

1.9000 

5.0739 

83.87 

.6255 

2.0000 

5.2259 

90.93 

.6172 

2.1000 

5.3780 

98.25 

.6095 

2.2000 

5.5300 

105.84 

.6021 

2.3000 

5.6820 

113.69 

.5952 

2.4000 

5.8340 

121.81 

.5886 

2.5000 

5.9860 

130.19 

.5823 

2.6000 

6.1380 

138.84 

.5764 

2.7000 

6.2901 

147.75 

.5707 

2.8000 

6.4421 

156.92 

.5653 

2.9000 

6.5941 

166.36 

.5602 

3.0000 

6.7461 

176.07 

.5553 
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